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Abstract 

Compact algebraic equations are derived, which connect the binding energy and the asymptotic nor- 
malization constant (ANC) of a subthreshold bound state with the effective-range expansion of the corre- 
sponding partial wave. These relations are established for positively-charged and neutral particles, using 
the analytic continuation of the scattering (S) matrix in the complex wave-number plane. Their accuracy 
is checked on simple local potential models for the 16 0+n, 16 0+p and 12 C+a nuclear systems, with exotic 
nuclei and nuclear astrophysics applications in mind. 
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There is a renewal of interest today for the quantum description of the low-energy scattering of 
two particles. This interest is mostly triggered by condensed-matter and ultracold-gases physics 
but similar studies are performed in other fields. Nuclear physicists in particular are studying low- 
energy collisions in the context of nuclear astrophysics and of exotic nuclei. A problem specific 
to nuclear physics is that positively-charged particles repel each other, which makes cross sections 
extremely small at low energy, and hence hard or impossible to measure. Theory thus plays an 
important role here but is made more complicated by the Coulomb interaction. 

Particular problems occur in the presence of weakly bound states: because of the wave nature 
of quantum phenomena, these play a role similar to resonances, leading to huge variations of cross 
sections. Such subthreshold bound states occur both in neutral cases (e.g., magnetic Feshbach 
resonances in atom-atom collisions or the "historical" deuteron bound state) and in charged cases. 
Famous examples are the lowest 2 + and \~ states of the 16 nucleus which, lying just below 
the 12 C+a threshold, are thought to strongly influence the 12 C(a,7) 16 capture cross section and 
hence the carbon to oxygen ratio in red giant stars yj, 120 . Essential quantities required for the 
theoretical description of this reaction are the binding energy of these states and their ANC, which 
characterizes the tail of their slowly-decreasing wave function [see Eq. (|T2|) below]. While the 
energy can be precisely measured, the ANC is still rather poorly known: it is not directly accessible 
experimentally and various indirect methods have been proposed to infer it. Among them, a high- 
precision measurement of elastic scattering is believed to give crucial information [3] but the ANC 
extraction from this measurement relies on a reaction-matrix analysis 040 , which is made rather 
delicate by the description of the non-resonant cross section pD. 

In the present work, a more fundamental approach is proposed to relate scattering properties 
to bound-state properties, in particular to the ANC, in the case of a weakly bound state. This 
approach is based on general S-matrix properties flg]: non-resonant scattering states are described 
with the help of the effective-range expansion while bound states are described in terms of poles 
of the S matrix in the complex plane. Combining these allows us to derive compact algebraic 
equations which prove particularly useful and promising, as illustrated by the examples below. In 
the following, we mostly concentrate on the charged case but also discuss the simpler neutral case, 
which leads to interesting comparisons. 

Let us consider two particles of charges Z\, Z 2 and of reduced mass fi. We denote the center- 
of-mass energy by E = h 2 k 2 /2fi and the dimensionless Sommerfeld parameter by 77 = l/a B k, 
where k is the relative wave number and a B = h 2 /ZiZ 2 e 2 n is the nuclear Bohr radius. The wave 
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number is allowed to be complex, the upper-half complex plane corresponding to the physical 
energy shee, g The scattering matrix S, for a partial wave , in the presence of hoth a Conionrb 
and a short range (e.g., nuclear) interaction is defined by [8Q 
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In these expressions, the first factor is the pure Coulomb scattering matrix, with ai(k) = arg r(Z + 
1 + irj) being the Coulomb phase shift. The second factor is due to the additional short-range 
interaction, with Si being the additional phase shift. 

In the following, we are interested in low physical energies (either positive or negative) and 
hence would like to expand the functions of interest around the origin E — 0. Since the S matrix 
defined above has a rather complicated analytical structure in the complex energy plane, due to 
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where two alternative expressions are given for both factors. In the last equation, we have used 
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Next, we define on the physical energy sheet H9I.I11H 

h(i] 2 ) = ip(ir]) — ln(i^) + (2i r q)^ 1 , 



where ip is the digamma function [12]. This allows constructing the effective-range function 
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which is holomorphic in the physical energy sheet, with a cut along part of the negative real axis 
BSD. This function is analytic and regular at the origin, which implies it can be Taylor expanded 
[see Eq. (fT6l) below]. Moreover, it is real for real energies. 

In the neutral case, things are much simpler: by taking i] = in Eqs. (0 and ©, the Coulomb 
factors simplify and the additional phase shift reduces to the usual phase shift. An effective-range 
function, with the same properties as in the charged case, can be defined as |7|] 
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Let us now connect the effective-range function with physical properties for both positive and 
negative energies. Forj^ositive energies, it is related to the scattering phase shift. In the charged 
case, Eq. © yields |9, 
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where the imaginary part has vanished because of the value of the digamma function for an imag- 
inary argument Jl?!] . In the neutral case, Eq. © implies \1 
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For negative energies, bound states correspond to poles of the S matrix on the positive imagi- 
nary k axis QVU- Let k = 1Kb be the location of such a pole, E b = —% 2 kI/2[i and r] b = —ijas^b be 
the corresponding binding energy and Sommerfeld parameter. Equation © shows that this bound 
state corresponds to a zero of Ff 1 and Eq. © implies that 111 ill 
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In the neutral case, Eq. © implies instead that 
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These equations relate the effective-range function to the bound-state energy. Though Eq. (ITTT) is 
frequently used in the literature, in particular to treat the I = deuteron bound state [7], Eq. (flOl) 
is less known. 

The ANC C& of this bound state is defined by the asymptotic behavior of its normalized radial 
wave function Rb(r), i.e., 

R b {r) ~ C b exp{-n b r)/r M+1 . (12) 
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In Ref. [| 14f | , this relation is shown to be valid for potentials decreasing faster than exp(— 2K b r). It 
thus holds, e.g., for all bound states of cut-off potentials, and for weakly bound states of potentials 
decreasing fast enough, a condition assumed in the present work. In other cases, in particular for 
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deeply bound states of exponentially decreasing potentials, Eq. (TT3l) can be violated [5]. Combin- 
ing Eqs. (00), ©, and CE3 implies that 111 ill 
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which provides a second relationship between the effective-range function and the bound- state 
properties. 

Equations (flOl) . (fTT|) . and (fl4l) are quite general; let us now particularize them to small binding 
and scattering energies by using the Taylor expansions of the various holomorphic functions in- 
troduced above for k 2 — > (i.e., for rf — > oo). The expansion of function wi directly ensues from 
its definition ©; function h expands as 
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and we write the first terms of function K[ as 

1 

ai + 2' 

The pole location conditions (flOl) and (TTT)) lead to different relations for the charged and neutral 
cases. Truncating Eq. (flOl) to the first order in E leads to 

1_ ( r, 1 

ai 

while truncating Eq. (fTTT) leads to 
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Both equations directly relate the scattering length ai to the binding energy, provided the effective 
range r\ is known. They imply that the scattering length is large when the binding energy is small, 
ri being in general rather independent of the binding energy. Note that, for I = 0, Eq. ([TBI is 
the well-known Sch winger relation [15]; for a very small binding energy, the second term then 
dominates, which implies that the scattering length becomes independent of the effective range. 
For I > 0, on the contrary, the first term dominates and Eq. (TT8l) becomes closer to Eq. (fT71) . 

Combining these equations with Eq. (fl4l) provides useful expressions for the ANC in terms of 
the scattering length. For the charged case, one gets 
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FIG. 1: ANC of a schematic 16 0+p bound state as a function of its binding energy, for the I = 0, 1 and 2 
partial waves and three different values of the potential diffuseness a. The two 1 F physical bound states 
are indicated by arrows. 



which generalizes Eq. (45) of Ref. |6D to I > 0. This equation is actually also valid for the neutral 
case with I > 0, where the fraction of the first member simplifies. For / = 0, one gets the very 
different result |7|] 
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These simple equations clearly illustrate the possibility of directly extracting the ANC from the 
binding energy and the elastic- scattering phase shifts. However, let us stress that their validity is 
limited to low energies, since they rely on a first-order effective-range expansion in energy. Higher- 
order expansions can be obtained [see, e.g., Eq. (f21~)) below] but lead to rather heavy expressions 
not reproduced here for brevity. Above a few terms, a purely numerical treatment is probably most 
appropriate. Let us also note that the precision on the ANC expected from these equations strongly 
depends on the precision available on the scattering length. In the presence of a subthreshold bound 
state, we have seen that the scattering length is very large. Since it appears in the effective-range 
expansion (fl~6l) with a negative power, it has a much smaller influence on this expansion than the 
following terms. Hence, it is typically difficult to determine ai from scattering phase shifts in such 
a case. Equations like (fTTT) and (fl~8l) then come as a useful tool: they allow one to calculate ai 
from the binding energy and the next terms of the expansion, which are easier to determine from 
scattering experiments. 
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FIG. 2: Same as Fig.Q] but for the 16 0+n system. 



Let us now check the accuracy of these results on typical nuclear physics examples. We first 
study the binding-energy dependence of ANCs for Woods-Saxon potentials of radius 3.023 fm, 
with and without Coulomb interactions. These potentials schematically represent the 16 0+p and 
16 0+n systems. For several values of the potential diffuseness, we vary the depth of these poten- 
tials, in order to scan a wide range of binding energies. For each binding energy, we calculate 
numerically the ANC. The results for the s, p and d waves are presented in Figs.[T|for the charged 
case and [2] for the neutral case. The noticeable simplicity of these curves up to about 1 MeV can 
be easily explained by Eqs. (fT71)- ((20l) : the ANCs have the expected energy dependence, showing 
the small sensitivity of the effective range r\ to the potential depth. On the contrary, its sensitivity 
to the potential diffuseness is larger and increases with the magnitude of the barrier due to the 
Coulomb and/or centrifugal potentials. For the remarkable s-wave neutral case [Eq. (|2Q|)1. the be- 
havior at small binding energy becomes independent of the potential, a situation also encountered 
for the deuteron bound state fl 1 5M . 

The 17 F and 17 nuclei both have an s and a d bound state, well represented by the above 
potentials. Their binding energies are indicated by arrows in Figs. [Hand |2l For the 17 F nucleus, 
the s state is weakly bound, which suggests a one-proton halo structure. Quantitatively, Eq. (fT9l) 
provides the ANC of this state with an accuracy better than 0. 1 %. Equation (fTTT) predicts ao with an 
error of 5%, while a generalization to the second order in energy reduces this error to 0.4%, which 
shows the fast convergence of our method. Now, for this system, a second-order effective-range 
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expansion turns out to be sufficient to fit experimental scattering data on a rather wide energy 
range (a few MeV) [16]. The present results thus show that both the scattering length and the 
bound-state ANC could be precisely deduced from these data. 

Let us now turn to a more important example from the physical point of view. The first 2 + 
excited state of 16 0, which lies 245 keV below the 12 C+a threshold, plays a key role in the cal- 
culation of the 12 C(a, 7) 16 radiative capture cross section jl, 2]. This state is known to have a 
d-wave 12 C+a structure and is reasonably well described by a 12 C+a potential model [5]. Here, we 
choose a nuclear Gaussian potential inspired by Ref. [17], which reads —112.3319 exp(— r 2 /2.8 2 ) 
MeV, where r is the internucleus distance in fm, and a Coulomb potential 12e 2 erf (r/2.5)/r, where 



erf is the error function [12]. This potential has a bound state at E b = —245.0 keV, corresponding 
to K b = 0.1876 fm" 1 with fr 2 /2/i = 6.964 MeV fm 2 and a B = 0.8060 fm. Numerically, one finds 
the ANC C b = 1.384 x 10 5 fm~ 1/2 . Let us now try to recover this value with Eq. (PT9| ). The numer- 



ical method of Ref. Q18D provides a 2 = 5.891 x 10 4 fm 5 , r 2 = 0.1580 fm" 3 , and P 2 = -65.96 fm 8 
for this potential. Equation (03 then leads to C b ps 1.23 x 10 5 fm~ 1/2 (11% error), while Eq. (flTT) 
provides a 2 ~ 4.76 x 10 4 fm 5 (19% error). A generalization of these equations to higher orders 
thus seems necessary in this case. Up to the third order in energy, the ANC formula reads 

icy . *M±MB (2d 
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By neglecting the last two terms of the denominator, one recovers Eq. (U9l) . Adding the second 
term leads to C b ~ 1.414 x 10 5 fm~ 1//2 (3% error). To determine the third term, we need the next 
coefficient of the effective-range expansion; for that, we have fitted the numerical phase shifts with 
Eqs. ([8]) and (TT6l) . which provides Q 2 = 0.1411(2) fm. Adding this term strikingly improves the 
fit: within a 1 % precision, the phase shifts are fitted up to 80 keV without it, and up to 700 keV 
with it. This improvement is very promising as experimental data, which are only available above 
750 keV |Q], could be fitted with a limited number of terms in the effective-range expansion. 
With the three terms, Eq. (|2TT) provides C b w 1.379 x 10 5 fm" 1 / 2 (0.4% error), a remarkable 
precision for a third-order effective-range expansion. 

In conclusion, starting from general S-matrix properties, we have established equations con- 
necting the effective-range function with the energy and ANC of bound states, for an arbitrary 
partial wave and for both charged and neutral systems. By using the effective-range series expan- 
sion, we have deduced from these relations compact equations [Eqs. (fl~7l) - (f2~Tj) . possibly extended 
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to higher orders] that relate the energy and ANC of a subthreshold bound-state to the series coef- 
ficients (scattering length, effective range. . . ). We have checked on test cases that these relations 
are precise and could be used to extract a subthreshold bound state ANC from scattering data in 
a model-independent way when precise experimental phase shifts are available. We plan to apply 
our method to the 12 C+a experimental phase shifts of Ref . [3] in a near future. 
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